SLAC-PUB-8824 
BABAR-PROC-01/27 
May, 2001 



MEASUREMENT OF CHARMLESS HADRONIC B DECAYS 
BRANCHING FRACTION AT BABAR. 

G.Cavoto 

on behalf of the BaBar Collaboration 
Universita di Roma La Sapienza, Dipartimento di Fisica and INFN, 

1-00185 Roma, Italy 
E-mail: gianluca.cavoto@romal.infn.it 



Abstract 

We present preliminary measurements of the branching fractions for charmless hadronic decays of 
B mesons into two-body final states with kaons, pions and a cj) resonance. The measurements are 
based on a data sample of approximately 23 million BB pairs collected by the BaBar detector at 
the PEP-II asymmetric B Factory at SLAC. 
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Charmless hadronic final states play an important role in the study of CP-violation. In the 
Standard Model, all CP-violating phenomena are a consequence of a single complex phase in the 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix El. Measurements of the rates and CP 
asymmetries for B decays into the charmless final states tttt and Kir can be used to constrain 
the angles a @ of the unitarity triangle. Decays containing a 4> meson are interesting since they 
are dominated by b — > s(d)ss peguin diagram, with potential benefits to estimates of direct CP- 
violation effects. They also allow an independent measurement of sin 2(3. 

We present preliminary measurements of the branching fractions for charmless hadronic decays 
of B mesons in the final states vr+vr", K+tt~, K + ir°, K°tt + , K°ir°, 4>K + , (f)K , and (892)1. 
The data sample used in these analyses consists of 22.57 x 10 6 BB pairs, collected at the PEP-II 
e + e~ collider (SLAC) with the BaBar detector! . 

Hadronic events are selected based on track multiplicity and event topology. We use only good 
quality tracks: tracks are identified as pions or kaons using the Cherenkov angle 8 C measured by 
a unique, internally reflecting Cherenkov ring imaging detector (DIRC). Candidate K° s mesons are 
reconstructed from pairs of oppositely charged tracks that form a well-measured vertex and have 
an invariant mass within 3.5 standard deviations (a) of the nominal K® massB. Candidate photons 
are defined as showers in the EMC that have the expected lateral shape and are not matched to a 
track. Candidate ir° mesons are formed from pairs of photons with an invariant mass within 3a of 
the nominal 7T° mass. The 7T° candidates are then kinematically fitted with their mass constrained 
to the nominal ir° mass. All tracks (except K® daughters) are required to have good-quality 9 C 
measurements that are inconsistent with the expected value for a proton. Electrons are rejected 
based on specific ionization (dE/dx) in the DCH system, shower shape in the EMC, and the ratio 
of shower energy to track momentum. 

Candidate B mesons are reconstructed in the various topologies: h + h'~, /i + 7r°, K®h + , K^tt , 
4>h + , (f>Kg, and cf>K* + where the symbols h and h' refer to ir or K. For eft candidates, both tracks 
must be identified as kaons whose invariant mass must lie within a ±30MeV/c 2 interval centered 
around (j) mass. The selection of K* + comprises K + tt° and K®ir + combinations within a Kit mass 
interval of ±150MeV/c 2 . A K* helicity angle cut effectively requires tt° momentum greater than 
0.35 GeV/c. 

Candidate B mesons are selected exploiting the kinematic constraints provided by the T(4S) 
initial state: we define a energy-substituted mass tties, where yfsjl is substituted for candidate's 
energy, and the difference AE between the B-candidate energy and yfs/2. For all modes the tties 
resolution is dominated by the beam energy spread and is approximately 2.5MeV/c 2 , while AE 
resolution is mode dependent and dominated by momentum resolution. Candidates are selected in 
the range 5.2 < toes < 5.3GeV/c 2 . Candidates are accepted, depending on the decay topology, in 
various AE ranges, which are restrictive enough to suppress backgrounds due to other types of B 
decays. 

The largest source of background is from random combinations of tracks and neutrals produced 
in the e + e~ — > qq continuum (where q = u, d, s or c). In the CM frame this background typically 
exhibits a two-jet structure. In contrast, the low momentum and pseudoscalar nature of B mesons 
in the decay T(4S) — * BB leads to a more spherically symmetric event. This topology difference is 
exploited using two event-shape quantities: the angle 6 S I between the sphericity axes, evaluated 
in the CM frame, of the B candidate and the remaining tracks and photons in the event. The 
distribution of the absolute value of cos# s is strongly peaked near 1 for continuum events and is 
approximately uniform for BB events. We require | cos# s | < 0.9 

The second quantity used in the analyses is a Fisher discriminant T which consists of a linear 



combination of several variables that distinguish signal from background 0. The experimental 
observables used in the definition of T are the scalar sum of the momenta of all tracks and photons 
(excluding the B candidate daughters) flowing into nine concentric cones centered on the thrust 
axis of the B candidate, in the CM frame. Each cone subtends an angle of 10° and is folded to 
combine the forward and backward intervals. 

The global detection efficiencies, which include the intermediate particle branching fractions, 
are listed in Table |l[ 

Table 1: Summary of results for detection efficiencies (e), numbers of fitted signal yields (Ns), statistical significances, 
and measured branching fractions (B). ( 90% confidence level upper limits). The efficiencies include the branching 
fractions for K° -> K% -> tt+tt", tt° -» 77, <j> -> K+ K~ . Equal branching fractions for T(4S) -» B°B° and B+ B~ 
are assumed. 



Mode 


e(%) 


N s 


Stat. Sig. (cr) 


£(HT e ) 




45 


41 ± 10 


4.7 


4.1 ± 1.0 ±0.7 


K+ir- 


45 


169 ± 17 


15.8 


16.7 ±1.6 ±1.3 


K+K- 


43 


n+7.8 


1.3 


< 2.5 


TT+TT 


32 


37 ± 14 


3.4 


< 9.6 


K+7T 


31 


75 ± 14 


8.0 


10.8ii;sj ± L0 


K°7T+ 


14 


oy -10 


9.8 


18.2±l;g ± 1.7 


K°K+ 


14 


0.0+2- 4 





< 2.5 


K°7T° 

(f>K + 


9.6 
18 


17-9±I;I 
31-4±g : S 


4.5 
10.5 


8.2±l^ ± 1.1 
7.7tii±0.8 




19 




0.6 


< 1.4 




6 




6.4 


8.1+y ±0.8 


(j)K* + 


5 




4.5 


9.7113 ± L7 


<pK* K + 


2.5 
2.4 




2.7 
3.6 


12.8+^ ±3.2 
8.0i|;§ ± 1.3 



Signal yields are determined from an unbinned maximum likelihood fit using mES; Ai2, JT, ^ 
mass, and X* + mass (where applicable). In each of the fits, the likelihood for a given candidate j is 
obtained by summing the product of event yield nk and probability Vk over all possible signal and 
background hypotheses k. The n& are determined by maximizing the extended likelihood function 
C: 

/ M \ N r M 

£ = exp -5>J JJ ^2n k P k (xj;a k ) . (1) 

V i=i / j=i U=i 

where Vk(xj]dk) is the probability for candidate j to belong to category (of M total categories), 
based on its characterizing variables xa and parameters a k that describe the expected distributions 
of these variables. The probabilities Vk(%j] <3fc) are evaluated as the product of probability density 
functions (PDFs) for each of the independent variables Xj, given the set of parameters d k - Monte 
Carlo simulated data is used to validate the assumption that the fit variables are uncorrelated. The 
exponential factor in C accounts for Poisson fluctuations in the total number of observed events N . 

The probabilities for each possible signal and background hypotheses are evaluated as the prod- 
uct of probability density functions (PDFs) for each of the independent variables. The parameters 
of 777. Eg, AE, and T PDFs are determined from data, and are cross-checked with Monte Carlo 



> 20 



BaBar 



H ii + + ++y+ + i M + + ii + + + ^Wv 

C i i i !_i ±1 i la*. 




20 



5.2 



4.++ , i + 1 1- 1 -t-f -H+ i 



5.225 



5.25 
m 



5.275 



5.3 



ES 



(GeV/c ) 



Figure 1: The tties distributions for candidates which are (from top to bottom) -k + tt -like, K + n -like, K + K -like 

The curves represent fits to the distributions 



simulation. In particular, a sample of -D*-tagged D° — > K^tt + decays is used to parametrize the 
9c distribution for pion and kaon tracks as a function of momentum. The results of the fits for 
the various topologies are summarized in Table |l[ We find evidence for the decay B + — ► tt + tt° 
and measure a branching fraction of B(B + — * tt + tt ) = (5.li£jj ± 0.8) x 10~ 6 . However, the signal 
significance is insufficient to claim observation. A 90% confidence level upper limit is computed 
for this mode, as well as for B° — ► K + K~ B + — ► K°K + , and -B + — > 0tt + in the following 
manner. The upper limit on the signal yield for mode k is given by the value of for which 

n° 

Jo k £max dnfc/ /Q^^Cmaxdnfc = 0.90, where £ max is the likelihood as a function of n&, maximized 
with respect to the remaining fit parameters. 

The result is then increased by the total systematic error. The detection efficiency is reduced 
by its systematic uncertainty in calculating the branching fraction upper limit. The statistical 
significance of a given channel is determined by fixing the yield to zero, repeating the fit, and 
recording the change in — 21n£. 

Figure [l] shows the distributions in m^s for events passing the selection criteria, as well as 
requirements on Fisher value, which are used to increase the relative fraction of signal events of a 
given type. Fits to these distributions are overlaid. 

Imperfect knowledge of the PDF shapes and of the detection efficiencies are the main sources 
of systematic uncertainties on the branching ratio measurements. Uncertainties in the PDF pa- 
rameterizations are estimated either by varying the PDF parameters within la of their measured 
uncertainties or by substituting alternative PDF from independent control samples, and recording 
the variations in the fit results. 

In summary, we have measured branching fractions for the rare charmless decays B° — ► ir + ir~ , 
B° -> K+tt-, B + -» K+tt , B+ -» K°tt + , B° -> K°tt°, (j)K + , 4>K°, and cf)K* + and set upper 
limits on B° -» K + K~, B+ -» tt+tt , B + -> K°K + , and 4>tt + . 
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